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Unstable separatrix solutions for the first and second Painleve transcendents are studied 
both numerically and analytically. For a fixed initial condition, say y( 0) = 0, there is a 
discrete set of initial slopes y'{ 0) = b n that give rise to separatrix solutions. Similarly, for 
a fixed initial slope, say y'{ 0) = 0, there is a discrete set of initial values y( 0) = c n that 
give rise to separatrix solutions. For Painleve I the large-n asymptotic behavior of b n is 
b n ~ Bin 3 ' 5 and that of c n is c n ~ C\n 2 / 5 , and for Painleve II the large-n asymptotic 
behavior of b n is b n ~ Bun 2 ^ 3 and that of c„ is c n ~ Cun 1 ^ 3 . The constants B\, C\, Bn, 
and Cn are first determined numerically. Then, they are found analytically and in closed 
form by reducing the nonlinear equations to the linear eigenvalue problems associated with 
the cubic and quartic 'PT-symmetric Hamiltonians H = |p 2 + 2 ix 3 and H = ^p 2 — \x 4 . 

PACS numbers: 02.30.Hq, 02.30.Mv, 02.60.Cb 


I. INTRODUCTION 

The famous Painleve transcendents are six nonlinear second-order differential equations whose 
key features are that their movable (spontaneous) singularities are poles (and not, for example, 
branch points or essential singularities). There is a vast literature on these remarkable differential 
equations m- These equations have arisen many times in mathematical physics; for a small 
sample, see Refs. IMS!. This paper considers the first and second Painleve transcendents, referred 
to here as P-I and P-II. The initial-value problem (IVP) for the P-I differential equation is 

y"(t) = 6[y(t)] 2 + t, y(0) = c, 2/'(0) = b (1) 

and the IVP for P-II (we have set an arbitrary additive constant to 0) is 

y"{t) = 2[2/(i)] 3 + ty(t), 2/(0) = c, y'(0) = b. (2) 

Many asymptotic studies of the Painleve transcendents have been published, but in this paper 
we present a simple numerical and asymptotic analysis that to our knowledge has not appeared 
in the literature. This analysis concerns the initial conditions that give rise to special unstable 
separatrix solutions of P-I and P-II. Our asymptotic analysis verifies the numerical results given in 
this paper for P-I and P-II as well as some preliminary numerical calculations that were presented 
in an earlier paper on nonlinear differential-equation eigenvalue problems m- 

The main idea, originally introduced in Ref. H3, is that a nonlinear differential equation may 
have a discrete set of critical initial conditions that give rise to unstable separatrix solutions. These 
discrete initial conditions can be thought of as eigenvalues and the separatrices that stem from these 
initial conditions can be viewed as the corresponding eigenfunctions. The objective in Ref. m was 
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to find the large-n (semiclassical) asymptotic behavior of the nth eigenvalue. The general analytical 
approach that was proposed was to simplify the nonlinear differential problem to a linear problem 
that could be used to determine the leading asymptotic behavior of the eigenvalues as n —> oo. 

A toy model was used in Ref. m to explain the concept of a nonlinear eigenvalue problem. 
This model makes use of the elementary first-order differential equation problem 

y'(t) = cosfvrf y(t)], y( 0) = a. (3) 

It was shown that the solutions to this initial-value problem pass through n maxima before van¬ 
ishing like 1/t as t —> oo. As the initial condition a = y{ 0) increases past special critical values 
a n , the number of maxima jumps from n to n + 1. At these critical values the solution y{t) to 
0 is an unstable separatrix curve in the following sense: At values of y(0) infinitesimally below 
a n the solution merges with a bundle of stable solutions all having n maxima and when y(0) is 
infinitesimally above a n the solution merges with a bundle of stable solutions all having n + 1 
maxima. The challenge is to determine the asymptotic behavior of the critical values a n for large 
n. (This generic problem is the analog of a semiclassical high-energy approximation in quantum 
mechanics.) To solve this problem it was shown that for large n, the nonlinear differential equation 
problem ([3]) reduces to a linear one-dimensional random-walk problem. The random-walk problem 
was solved exactly, and it was shown analytically that 

a n ~ 2 5 ' /6 y / n (n —> oo). (4) 

Kerr subsequently found an alternative solution to this asymptotics problem and verified @ m- 
The nonlinear eigenvalue problem described above is similar in many respects to the linear 
eigenvalue problem for the time-independent Schrodinger equation. For a potential V(x) that rises 
as x —> Too, the eigenfunctions of the Schrodinger eigenvalue problem 

— ip"(x) + V(x)i/j(x) = Eip(x), '0(d=oo) = 0, (5) 

are unstable with respect to small changes in the eigenvalue E; that is, if E is increased or decreased 
slightly, ^(x) abruptly ceases to obey the boundary conditions [and thus is not normalizable (square 
integrable)]. Furthermore, like the eigenfunctions (separatrix curves) of ([3]), the eigenfunction 
ij} n {x) corresponding to the nth eigenvalue has n oscillations in the classically allowed region before 
decreasing monotonically to 0 in the classically forbidden region. 

This paper considers four eigenvalue problems. First, for P-I we find the large-n behavior of the 
positive eigenvalues b n for the initial condition y( 0) = 0, y'{ 0) = b n and also the large-n behavior 
of the negative eigenvalues c n for the initial condition y(0) = c n , y'{ 0) = 0. We show that 

b n ~ Bin 3 / 5 and c n ~ C\n 2 ^. 


Second, for P-II we show that for large n the asymptotic behaviors of b n and c n are given by 

b n ~ Bun 2 / 3 and c n ~ Cun 1 / 3 . 

We determine the constants B\, C\, Bu, and C\\ both numerically and analytically. 

This paper is organized as follows. In Sec. [TT] we obtain the constants B\ and C\ by using 
numerical techniques and in Sec. |III| we do so analytically by reducing the large-eigenvalue problem 
to the linear time-independent Schrodinger equation for the cubic "PT-symmetric Hamiltonian 
Next, we study the eigenvalue problem for the second Painleve transcendent. 


H = \p 2 + ix 3 


In Sec. IV we present a numerical determination of the large-n behavior of the eigenvalues and 
in Sec. [V] we verify the numerical results in Sec. IV by using asymptotic analysis to reduce the 
nonlinear large-eigenvalue problem for P-II to the linear Schrodinger equation for the quartic VT- 
symmetric Hamiltonian H = ^ p 2 — \x A . In Sec. VI we make some brief concluding remarks. 
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II. NUMERICAL ANALYSIS OF THE FIRST PAINLEVE TRANSCENDENT 

In Ref. [17] there is a brief numerical study of the initial-value problem for the first Painleve 
transcendent ([Tj) . It is easy to see that there are two possible asymptotic behaviors as t —> — oo; the 
solutions to the P-I equation can approach either +y/—t /6 or —yj—t/ 6 . An elementary asymptotic 
analysis shows that if the solution y(t) approaches —y/—t/ 6 , the solution oscillates stably about 
this curve with slowly decreasing amplitude m- However, while the curve +yf—t/Q is a possible 
asymptotic behavior, this behavior is unstable and nearby solutions tend to veer away from it. 
We define the eigenfunction solutions to the first Painleve transcendent as those solutions that 
do approach +y/—t /6 as i -> —oo. These separatrix solutions resemble the eigenfunctions of 
conventional quantum mechanics in that they exhibit n oscillations before settling down to this 
asymptotic behavior. However, because the P-I equation is nonlinear, these oscillations are violent; 
the nth eigenfunction passes through [n/ 2 ] double poles where it blows up before it smoothly 
approaches the curve +y/—t/6. (The symbol [n/2] means greatest integer in n/2.) 

One can specify two different kinds of eigenvalue problems for P-I, each of which is fundamentally 
related to the instability of the asymptotic behavior +y/—t/6. One can (i) fix the initial value y( 0) 
and look for (discrete) values of the initial slopes y'( 0 ) = b that give rise to solutions approaching 
+y/—t/6, or else (ii) one can fix the initial slope y'(Q) and look for the (discrete) initial values of 
y{0 ) = c that give rise to solutions approaching +yj—t/d>. 


A. Initial-slope eigenvalues for Painleve I 

Let us examine the numerical solutions to the initial-value problem for the P-I equation ([Tj) 
for t < 0. To find these solutions we use Runge-Kutta to integrate down the negative-real axis. 
When we approach a double pole and the solution becomes large and positive, we estimate the 
location of the pole and integrate along a semicircle in the complex-t plane around the pole. We 
then continue integrating down the negative-real axis. We choose the fixed initial value y(0) = 0 
and allow the initial slope y'{ 0) = b to have increasingly positive values. (We only present results 
for positive initial slope; the behavior for negative initial slope is analogous and describing it would 
be repetitive.) Our numerical analysis shows that the particular choice of y( 0) is not crucial; for 
any fixed y( 0 ) the large-n asymptotic behavior of the initial-slope eigenvalues b n is the same. 

We find that above the critical value bi = 1.851854034 (the first eigenvalue) there is a continuous 
interval of b for which y(t) first has a minimum and then has an infinite sequence of double poles 
(see Fig. [TJ left panel). However, if b increases past the next critical value &2 = 3.004031103 (the 
second eigenvalue), the character of the solutions changes abruptly and y(t) oscillates stably about 
-yJ-t/G (Fig- 0 right panel). When b exceeds the critical value 63 = 3.905175320 (the third 
eigenvalue), the solutions again exhibit an infinite sequence of poles (Fig. [2j left panel). When b 
increases past the fourth critical value 64 = 4.683412410 (fourth eigenvalue), the solutions once 
again oscillate stably about 6 (Fig. [TJ right panel). Our numerical analysis indicates that 

there is an infinite sequence of critical points (eigenvalues) at which the P-I solutions alternate 
between infinite sequences of double poles and stable oscillation about — y/—t/ 6 . 

The solutions that arise when j/(0) is at an eigenvalue have a completely different (and unstable) 
character from those in Figs.jljan d! These special solutions pass through a finite number of double 
poles (analogous to the oscillatory behavior of quantum-mechanical bound-state eigenfunctions in 
the classically allowed region of a potential well) and then undergo a turning-point-like transition 
in which the poles cease and y(t) exponentially approaches the limiting curve +\J—t/ 6 . The 
solutions arising from the first and second critical points b\ and 62 are shown in Fig. [3j those 
arising from the third and fourth critical points 63 and 64 are shown in Fig. |4j and those arising 
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FIG. 1: Typical behavior of solutions to the first Painleve transcendent y(t) for the initial conditions y( 0) = 0 
and b = y'( 0). In the left panel b = 2.504031103, which lies between the eigenvalues &i = 1.851854034 and 
62 = 3.004031103. In the right panel b = 3.504031103, which lies between the eigenvalues 62 = 3.004031103 
and 63 = 3.905175320. The dashed curves are y = ± \J—t/ 6 . In the left panel the solution y(t) has an 
infinite sequence of double poles and in the right panel the solution oscillates stably about — \Jt/ 6. 



FIG. 2: Solutions to the P-I equation (jTJ) for y( 0) = 0 and b = y'( 0). Left panel: b = 4.583412410, which 
lies between the eigenvalues 63 = 3.905175320 and 64 = 4.6834124103. Right panel: b = 4.783412410, which 
lies between the eigenvalues 64 = 4.683412410 and 65 = 5.383086722. 


from the tenth and eleventh critical points 610 and bn are shown in Fig. [5} The critical points are 
analogous to eigenvalues because they give rise to unstable separatrix solutions; if y'{ 0 ) changes 
by an infinitesimal amount above or below a critical value, the character of the solutions changes 
abruptly and the solutions exhibit the two possible generic behaviors shown in Figs. [T] and [2] 

In Ref. [17] a numerical asymptotic study of the critical values b n for n 3> 1 was performed by 
using Richardson extrapolation [20]. [In Ref. [T7] the initial value was chosen to be 2 /( 0 ) = 1 rather 
than 2 /( 0 ) = 0 as in the current paper. However, as emphasized above, if 2 /( 0 ) is held fixed, we find 
that the large-n asymptotic behavior of the initial slope b n is insensitive to the choice of 2/(0)-] It 
was found in Ref. m that for large n, the nth critical value had the asymptotic behavior 

24 ( 0 ) = b n ~ .Bin 3 / 5 (n 00 ). ( 6 ) 

In Ref. [T7] the constant B\ was determined numerically to an accuracy of about four or five decimal 
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FIG. 3: First two separatrix solutions (eigenfunctions) of Painleve I with initial condition y{ 0) = 0. Left 
panel: y'( 0) = b\ = 1.851854034; right panel: y'{ 0) = b 2 = 3.004031103. The dashed curves are y = 

±V^i/6. 



FIG. 4: Third and fourth eigenfunctions of Painleve I with initial condition y{ 0) = 0. Left panel: y'{ 0) = 
b 3 = 3.905175320; right panel: y'{ 0) = 64 = 4.683412410. 


places. However, we have now performed a more accurate numerical determination of the constant 
B\ by applying fifth-order Richardson extrapolation to the first eleven eigenvalues, and we have 
found the value of B\ accurate to one part in nine decimal places: 

Bi = 2.09214674. (7) 

On the basis of our numerical analysis, we can say with confidence that the underlined digit lies 
in the range from 3 to 5, so our determination of B\ is accurate to one part in 2 x 10 8 . 


B. Initial-value eigenvalues for Painleve I 

If we hold the initial slope fixed at y'{ 0) = 0 and allow the initial value y( 0) = c to become 
increasingly negative, we find that there is a sequence of negative eigenvalues c n for which the solu¬ 
tions behave like the eigenfunction separatrix solutions in Figs. [3]^5} The first four eigenfunctions 
are plotted in Figs. [ 6 ] and [7J 
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FIG. 5: Tenth and eleventh eigenfunctions of Painleve I with initial condition y(0) = 0. Left panel: 2/(0) = 
bio = 8.244932302; right panel: y'{ 0) = bn = 8.738330156. Note that as n increases, the eigenfunctions pass 
through more and more double poles before exhibiting a turning-point-like transition and approaching the 
limiting curve +y/—t/6 exponentially rapidly. This behavior is analogous to that of the eigenfunctions of 
a time-independent Schrodinger equation for a particle in a potential well; the higher-energy eigenfunctions 
exhibit more and more oscillations in the classically allowed region before entering the classically forbidden 
region, where they decay exponentially. 




FIG. 6: First two separatrix solutions (eigenfunctions) of Painleve I with fixed initial slope 2/(0) = 0. 
Left panel: 2/(0) = Ci = —0.7401954236; right panel: 2/(0) = c-i = —1.206703845. The dashed curves are 
2 / = 


Applying fourth-order Richardson extrapolation to the first 15 eigenvalues, we find that for 
large n the sequence of initial-value eigenvalues c n is asymptotic to Cin 2 / 5 , where the numerical 
value of the constant C\ is 


Ci = -1.0304844. (8) 

We are confident that the final digit is accurate to an error of ±1 and thus Ci is determined to an 
accuracy of one part in 10 7 . 
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3 >- 0 



FIG. 7: Third and fourth eigenfunctions of Painleve I with initial slope y'( 0) = 0. Left panel: y( 0) = C 3 
— 1.484375587; right panel: y( 0) = c 4 = -1.69951765. 


III. ASYMPTOTIC CALCULATION OF Bj AND Ci 


In this section we present an analytic calculation of B\ and C\ in ([7]) and Q. To begin, we 
multiply the P-I differential equation in (jTj) by y'{t) and integrate from t = 0 to t = x. We get 

H ee i[y'(*)] 2 - 2[y( X )f = Ifo'fO)] 2 - 2[y(0)] 3 + I(x), (9) 

where I(x) = f ( J dtty'{t). Note that the path of integration is the same as that used to calculate 
y(t) numerically in Sec. [TIJ it follows the negative-real axis until is gets near a pole, at which point 
it makes a semicircular detour in the complex-f plane to avoid the pole. 

If we evaluate I{x) for large |x| in the classically allowed region (just before the poles abruptly 
cease at the turning point), we find that as n -» oo, I(x) fluctuates and becomes small compared 
with H. This is not surprising because I(x ) receives many positive and negative contributions from 
the poles. [In fact, by calculating I(x) as x —> —oo, we can see a clear signal of an eigenvalue; as 
y'{0 ) = b passes an eigenvalue, I(x) goes from having positive to negative (or negative to positive) 
fluctuations but at an eigenvalue I(x) is smooth and not fluctuating.] Thus, for large n we treat 
the fluctuating quantity I(x ) as small, and if we do so we can interpret H as a time-independent 
quantum-mechanical Hamiltonian. [The isomonodromic properties of H when I (x) is not neglected 
were studied in Ref. [6].] 

We conclude that the large-n (semiclassical) behavior of the eigenvalues [that is, the initial 
conditions in ([!])] can be determined by solving the linear quantum-mechanical eigenvalue problem 
Hip = Eif>, where H = ^ p 2 — 2x 3 . To find these eigenvalues we rotate H into the complex plane 
[21] and obtain the well-studied PT-symmetric Hamiltonian [22] 

H = r^p 2 + 2 ix 3 . (10) 


The large eigenvalues of this Hamiltonian can be found by using the complex WKB techniques 
discussed in detail in Ref. |22j . For the general class of PT-symmetric Hamiltonians H = \p 2 + 
gx 2 {ix) e (e > 0), the WKB approximation to the nth eigenvalue (n 1) is given by 


En ~ s) 2/(4+e ) 


r (i + Fh) \^ n 

sm (jfs) T (l + 


(2e+4)/( e +4) 


( 11 ) 
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Thus, for H in (10) we take g = 2 and e = 1 and obtain the asymptotic behavior 

2 


E„ 


g/5 

VSnT (^) n/T (|) (n—»oo). 


( 12 ) 


Since if in (10) is time independent, we can evaluate H in (|9j) for fixed y(0) and large y'( 0) = b n 

(13) 


and obtain the result that 

b n ~ \pZEn = 5m 3 / 5 (n —» 00 ), 

which verifies <§• We then read off the analytic value of the constant B\\ 

Bi = 2 


V^r($) /r (i)l 3/5 


(14) 


which agrees with the numerical result in ([ 7 ]). Also, if we take the initial slope y'(0) to vanish and 
take the initial condition y(0) = c n to be large, we obtain an analytic expression for Ci, 


Ci = - 

which agrees with the numerical result in 


^r($)/r($)l 2/5 


(15) 


IV. NUMERICAL ANALYSIS OF THE SECOND PAINLEVE TRANSCENDENT 


To understand the behavior of solutions to the initial-value problem in ([ 2 ]) for Painleve II, we 
follow the procedure used in Sec. [IT] to study P-I. An elementary asymptotic analysis shows that as 
t —> — 00 , there are three possible asymptotic behaviors for solutions y(t). First, y(t) can oscillate 
stably about the negative axis. Second, y(t) can approach the curves ±y/—t/ 2; however, both of 
these asymptotic behaviors are unstable. 

If we numerically integrate ([ 2 ]), we observe that when t becomes large and negative, a typical 
solution to the P-II initial-value problem either oscillates about the negative axis or passes through 
an infinite sequence of simple poles. However, it is also possible to find special eigenfunction so¬ 
lutions that pass through only a finite number of poles and then approach either the positive or 
the negative branches of the square-root curves. These eigenfunctions obey the boundary condi¬ 
tions y{ 0) = 0 and y'{ 0) = ±b. [Note that P-II is symmetric under y — > —y, so there are two 
sets of eigenfunctions, one for each sign of y'{ 0).] We study these eigenfunctions numerically in 
Subsec. IV A The P-II equation is particularly interesting because as t — > + 00 , the behavior y — > 0 
becomes unstable. Thus, it is possible to have new kinds of eigenfunctions for positive t as well. We 
seek eigenfunctions that satisfy y'{ 0) = 0 and y{ 0) = c and examine the positive-c eigenfunctions 


numerically in Subsec. IV B 


A. Initial-slope eigenvalues for Painleve I 

Similar to what we found in Sec. [nj if we choose y( 0) = 0, there are critical values y'{ 0) = b n at 
which the solutions y(t) change their character. In Figs. [ 8 ] and [9] we plot the solutions to the P-II 
equation for the initial condition y( 0) = 0 and y'( 0) = b for b\ < b < 62 , 62 < b < 63, 63 < b < 64, 
and 64 < b < 65. Note that in these figures the character of the solution alternates between having 
an infinite sequence of simple poles and oscillating stably about y(t ) = 0. However, when y'{ 0) = b 
is at a critical value (eigenvalue) b n , the solution y(t) passes through a finite number [n/ 2 ] of simple 
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FIG. 8 : Typical behavior of solutions to the second Painleve transcendent for the initial conditions y{ 0) = 0 
and b = y'{ 0). In the left panel b = 1.028605106, which lies between the eigenvalues &i = 0.5950825526 and 
&2 = 1.528605106. In the right panel b = 2.028605106, which lies between the eigenvalues 62 = 1.528605106 
and 63 = 2.155132869. In the left panel the solution y(t) has an infinite sequence of simple poles and in the 
right panel the solution oscillates stably about — y/t/6. The dashed curves are the functions ± \J—tj 2. 




FIG. 9: Solutions to the P-II equation Q for y( 0) = 0 and b = y'( 0). Left panel: b = 2.600745985, which 
lies between the eigenvalues b 3 = 2.155132869 and 64 = 2.700745985. Right panel: b = 2.800745985, which 
lies between the eigenvalues 64 = 2.700745985 and 65 = 3.195127590. 


poles and then approaches either + y/—t/ 2 or 2. These eigenfunctions (separatrices) are 

and 12 for n = (1,2), (3,4), and (20,21). 

and [12] alternate between approaching the upper- 


plotted in Figs. 10, 11 


Note that the eigenfunctions in Figs. [lOj [TlJ 
unstable branch +y/—t/ 2 or the lower-unstable branch —\J—t/ 2 , and thus there are actually two 
sequences of eigenvalues, one for even n and one for odd n. Using Richardson extrapolation, we 
find that the sequences of eigenvalues & 2 n and & 2 n+i have the same asymptotic behavior 


b 2 n ~ b 2 n +1 ~ Bun 2 ^ 3 (n —> 00 ). 


(16) 


Our numerical calculations give 


Bn = 1.8624128. 


(17) 
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FIG. 10: First two separatrix solutions (eigenfunctions) of Painleve II with initial condition y( 0) = 0. Left 
panel: y'( 0) = b\ = 0.5950825526; right panel: y'( 0) = b 2 = 1.528605106. The dashed curves are ±\j~t/2. 



FIG. 11: Third and fourth eigenfunctions of Painleve II with initial condition y(0) = 0. Left panel: y'( 0) = 
b 3 = 2.155132869; right panel: y’{ 0) = 64 = 2.700745985. 


The numerical data for P-II are slightly more noisy than those for P-I, and fourth-order Richardson 
extrapolation only gives the underlined eighth digit as 8 ± 2 . 


B. Initial-value eigenvalues for Painleve II 


Next, we plot the positive-t solutions to P-II for vanishing initial slope and positive initial 
condition for t > 0. As t —> 00 , the nth eigenfunction passes through n simple poles before it 
approaches zero monotonically. In Figs. 13 14, and[T5]we plot the six eigenfunctions corresponding 


to n = (1,2), (3,4), and (13,14). (Because of the symmetry of P-II, for every positive eigenvalue 
there is a corresponding negative eigenvalue. We do not plot the negative-eigenvalue solutions.) 
Using fourth-order Richardson we determine that for large n, c n ~ Cun 1 / 3 , where 


Cu = 1.21581165. 


(18) 


The last digit 5 has an uncertainty of ±1. 
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FIG. 12: The twentieth and twenty-first eigenfunctions of Painleve II with initial condition y( 0) = 0. Left 
panel: y'( 0) = 620 = 8.499476190; right panel: y'( 0) = 621 = 8.787666814. 



FIG. 13: First two separatrix solutions (eigenfunctions) of Painleve II with fixed initial slope y'(0) = 0. Left 
panel: y{ 0) = C\ = 1.222873339; right panel: y( 0) = c 2 = 1.533883935. 


V. ASYMPTOTIC CALCULATION OF B n AND C n 


To obtain analytic expressions for Bn in ( flT| ) and Cn in ( [L8] ) , we follow the same procedure 
as in Sec. Ill for P-I. We multiply the P-II differential equation in ([ 2 ]) by y'{t) and integrate from 
t = 0 to t = x, where x is in the turning-point region which the simple poles stop. The result is 


H ee \[y\x)f - \[y{x)f = \[y\ 0)] 2 - ^(O)] 4 + /(*), 


(19) 


where I(x) = f Q l dtty{t)y'(t). The path of integration is the same as that used to calculate P-II 
numerically in Sec. |IV[ it follows the negative-real axis until it gets near a simple pole, at which 


point it makes a semicircular detour in the complex-! plane to avoid the pole. Again, as in Sec. Ill 


we argue that along this path the integrand of I(x) is oscillatory and because of cancellations we 
may neglect I(x) when n is large. 

We treat H as the VT -symmetric quantum-mechanical Hamiltonian 


H = 


( 20 ) 
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FIG. 14: Third and fourth eigenfunctions of Painleve II with initial slope 2/(0) = 0. Left panel: y{ 0) = C3 = 
1.754537281; right panel: 2/(0) = c 4 = 1.93061783. 



FIG. 15: Thirteenth and fourteenth separatrix solutions (eigenfunctions) of Painleve II with fixed initial 
slope 2 /( 0 ) = 0. Left panel: 2/(0) = c-\ = 2.858869051; right panel: 2/(0) = C 2 = 2.9303576515. 


and we use ( |11[ ) with g = 1/2 and e = 2 to obtain the formula 

. 1 4/3 


E r 


1 

~ 2 


3nv/2^r (I) /r (I) 


for the large eigenvalues of H. Finally, we calculate the eigenvalues b n by using 

1 ] 2 /3 


\f^En ~ |3nv&(|)/rQ) 
This result allows us to identify the value of B\\ in (JT7|) as 


(n -» 00 ). 


B u = 


3v^fr(|) /r(i) 


2/3 


( 21 ) 


( 22 ) 


(23) 


This result agrees with the numerical determination in ( |17[ ). 

To calculate Cjj we observe from Figs. 13 15 that the initial value y{ 0) is positive. However, if 
we neglext I{x) and assume a vanishing initial slope, we see that the right side of (19) negative. 















































































13 


Thus, as we did for the cubic Hamiltonian \p 2 — 2x 3 , we perform a complex rotation of the coupling 
constant to convert the quartic Hamiltonian to the form 


H = \f + \x\ 

(24) 

This is the conventional Hermitian quartic-anharmonic-oscillator Hamiltonian, and does not belong 
to the class of PT-symmetric Hamiltonians H = ^ p 2 + gx 2 (ix) e . A WKB calculation gives the 
large-eigenvalue approximation 

E n ~ [3ny/irT (§) /V (|)] 4/3 (n -> oo). 

(25) 

Thus, we read off the value of Cjj: 


c H = [3^r(f)/r(|)] 1/3 , 

(26) 


which agrees exactly with the numerical result in (|18[). 


VI. BRIEF CONCLUDING REMARKS 

In this paper we have shown that the first two Painleve equations, P-I and P-II, exhibit in¬ 
stabilities that are associated with separatrix solutions. The initial conditions that give rise to 
these separatrix solutions are eigenvalues. We have calculated the semiclassical (large-eigenvalue) 
behavior of the eigenvalues in two ways, first by using numerical techniques and then by using 
asymptotic methods to reduce the initial-value problems for the nonlinear P-I and P-II equations 
to linear eigenvalue problems associated with the time-independent Schrodinger equation. The 
agreement between these two approaches is exact. 

The obvious continuation of this work is to examine the next four Painleve equations, P-III 
- P-VI, to see if there are instabilities, separatrices, and eigenvalues for these equations as well. 
However, the techniques we have applied here may also be useful for other nonlinear differential 
equations such as the Thomas-Fermi equation y"(x) = [y(x)j 3 ^ 2 /y/x, which is posed as a boundary- 
value problem satisfying the boundary conditions y{ 0) = 1 and y( oo) = 0. The solution to this 
problem is unstable with respect to small changes in the initial data; if the initial slope y'{0 ) is 
varied by a small amount, the solution develops a spontaneous singularity at some positive value 
a. A leading-order local analysis suggests that this singularity is a fourth-order pole of the form 
400(x — a) -4 . However, this singularity is not a pole. Indeed, a higher-order local analysis indicates 
that there is a logarithmic-branch-point singularity at x = a as well and thus the solutions to the 
Thomas-Fermi equation live on multisheeted Riemann surfaces. It would be interesting to see if 
our work on nonlinear eigenvalue problems extends beyond meromorphic functions. 
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